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This paper collects miscellaneous results about the group §U(1, 1) that are helpful in applications 
in quantum optics. Moreover, we derive two new results, the first is about the approximability of 
SU(1, 1) elements by a finite set of elementary gates, and the second is about the regularization of 
group identities for tomographic purposes. 

I. INTRODUCTION 

In the last decades many achievements in quantum optics came from nonlinear effects in crystals (for a review on 
the topic see 0]). Nonlinear crystals allowed to produce both single mode squeezed states, which carry attenuated 
quadrature noise and constitute good carriers for classical information ||, [| |j, and two- modes squeezed states, such as 
the twin beam, which is a prototype for harmonic oscillator entangled states and are useful in many applications, such 
as continuous variables teleportation ||. On mathematical grounds, the action of nonlinear crystals can be described 
by parametric unitary transformations in which the pump mode is considered as a classical field and its creation and 
annihilation operators are substituted by the complex amplitude. The effective Hamiltonian allows the parametric 
down conversion, which is the process by which a photon with high frequency is annihilated and two photons with 
lower frequencies are created. This process gives rise to time evolution that can be described through unitaries in the 
Schwinger representation of the group SU(1, 1), namely exponentials of linear combinations of the three generators 

K + = cJb\ K_=ab, K z = -(a^a + tfb + l), (1) 

where a and b are the annihilation operators for the two modes. The degenerate parametric down conversion happens 
when the two created photons are in the same mode with a frequency which is half of the annihilated photon frequency, 
and this particular case giving rise to single mode squeezing corresponds to a — b with the three generators 

A' + = i(at)2, K_ = \a\ tf, = i(ata+l/2). (2) 

Squeezed states and twin beams are nowadays widely used in experimental quantum optics, and it is clear that 
the ability of manipulating radiation modes by unitaries of the group SU(1, 1) is crucial. In this paper we consider 
some general aspects of the group SU(1, 1) that can be exploited on the physical ground in order to approximately 
simulate any SU(1, 1) transformation by a finite set of elementary gates, namely unitary transformations which can 
be applied in a given succession in order to approach a target unitary in the representation of SU(1, 1). This is very 
useful in a situation in which an experimenter needs a flexible setup which allows to simulate within some accuracy 
any possible gate. A similar situation holds for qubits, where a very powerful theorem due to Solovay and Kitaev 
states that any gate can be efficiently approximated by a finite set of elementary gates. In the case of harmonic 
oscillators, however, the theorem still lacks an important part, which states that the amount of elementary gates 
needed in order to approximate any gate grows logarithmically with the accuracy. This fact is due to the dimension 
of the Hilbert space, and some intermediate result toward the analog of the qubit Solovay-Kitaev theorem in the case 
of harmonic oscillators can be derived with the reasonable assumption that the states of interest on which the gates 
have to be applied have finite average energy and finite variance of the energy distribution. In the paper we will also 
discuss severe limitations that forbid to find a power law which is independent of the group element that one wants 
to approximate. 

Besides the problem of approximation of squeezed states we can consider the problem of classifying and analyzing the 
performances of covariant measurements and tomographic measurements. The first ones are an idealization of physical 
measurements which turns out to be interesting because they saturate bounds on precision for the estimation of 
squeezing parameters, thus providing an absolute standard for rating of actual detectors. As regards the tomographic 
measurements, their statistics allows to completely determine the state of radiation modes — up to statistical errors. 
In Ref. H a particular tomographic measurement has been proposed for states with even or odd parity, based on 
properties of the Schwinger representation of SU(1, 1). In this paper we will discuss the possibility of deriving similar 
tomographic identities from group integrals. Moreover, an interesting mechanism because of which the "natural" 
group integral does not converge for physically interesting representations, and a sort of regularization is needed 
is shown. This analysis provides a whole range of tomographic POVMs corresponding to different regularizations, 
which can be studied in order to optimize the performances of SU(1, 1) tomography. The technique is general and 
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can be applied to many tomographic measurements originated from other groups. The core of the regularization 
technique consists in modifying the invariant (Haar) measure on the group manifold, and this modification gives rise 
to a generalization of the Duflo-Moore |t], || operator which is typical in groups which are not unimodular, namely 
for which the invariant Haar measure does not exist. This fact implies some complication in the data processing with 
respect to the usual homodyne tomography, but on the other hand allows to optimize the group measure in order to 
minimize the statistical errors. 

In Section |n| we discuss some general aspects of the group SU(1, 1), considering its defining representation. The 



results derived there will be exploited in subsequent sections. In Section III we prove the existence of a set of three 
elementary gates, and discuss the possibility to use them for approximation of target group elements under reasonable 
assumptions on the physical states. We also discuss the impossibility of having the exact analog of the Solovay-Kitaev 
theorem for the quantum optical representations of SU(1, 1). In Section [V] we show that the physical representations 
of SU(1, 1) are not square summable, and we show how one can modify the group theoretical identities for group 



integrals in order to obtain converging integrals which are useful for group tomography. In Section VI we close the 
paper with a summary of the contents and concluding remarks. 

II. GENERAL ASPECTS OF THE GROUP SU(1, 1) 

SU(1, 1) is the group of complex 2x2 matrices M with unit determinant that satisfy the relation 

M f PM = P , (3) 

where 

p =(;»)■ «> 

This relation implies that the elements of SU(1, 1) preserve the Hermitian form ui(vi, Ua) = v[Pv2 for arbitrary column 
vectors Vi € C 2 . 

From the above definition it is simple to show that any matrix M £ SU(1, 1) has the form 

for a, f3 complex numbers such that \a\ 2 — |/3| 2 = 1. Notice that the columns M\,M2 of M are orthogonal and 
normalized with respect to the form u, namely uj(Mi, M2) = 0, u>{Mx,M{) = 1, and u)(M2, M2) = — 1. By writing 
a = t + iz and f3 = x + iy, we obtain 

M = tl + iza z + xa x + y<r y , t 2 + z 2 - x 2 - y 2 = 1 , (6) 

1 and cr x ,<j y , a z being the identity and the three Pauli matrices, respectively. In other words, the elements of SU(1, 1) 
are parametrized by points of an hyperboloid in R 4 . This makes SU(1, 1) a Lie group, namely a group which is also 
a differentiable manifold. The above parametrization clearly exhibits three relevant facts: i) a group element is in 
one-to-one correspondence with three real parameters (x, y and z, for example), namely the group manifold is three 
dimensional ii) the group SU(1, 1) is not compact, and in) it is not simply connected [JL2|. 

Given a parametrization M(r), where the element M(r) £ SU(1, 1) is specified by the triple r £ R 3 , the matrix 
multiplication induces a composition law in the parameter space: (r, sjuros, where ro s is defined by the relation: 
M{r o s) = M(r)M(s). In particular, if r = (x, y, z), with x, y, z as in Eq. (^), we can define the invariant measure 

d(i(r) = 1 ==dxdydz . (7) 

y/ 1 + x z + y z — z A 

Invariance of the measure means that the action of the group does not change the volume of regions in the parameter 
space, namely, for any r,s, d/i(r o s) = d/i(so f) = dfi(r). The expression (0) of the invariant measure dfi(x, y, z) 
is particularly useful, since it allows to obtain the invariant measure in any parametrization of the group, just by 
performing a change of variables. For example, a useful alternative parametrization of a group element M £ SU(1, 1) 
is given by 

coshtf sinh6> 

e 



M (^)=[ diffe* eosh^-^ ) • ^ 
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for 6 £ [0, +oo), cj> £ [0, 2ir), ip £ [0, 2ir). The change of parametrization from (||) to (||) corresponds to the change of 
variables x = sinh cos tp ,y — sinh 8 sin ip , z — cosh 9 sin <f>. Performing the change of variables in Eq. (|t]) we obtain 
the expression of the invariant measure in the parametrization M — M{6, 4>, ip), namely 

di/(0, <p, i>) = sinh 6 cosh 6 d<9d<W . (9) 



A. The Lie algebra su(l, 1) 



Since SU(1, 1) is a real three-dimensional manifold, its Lie algebra su(l, 1) — the tangent space in the identity — is 
a three-dimensional vector space. As usual, a basis of the Lie algebra is obtained by differentiating curves passing 
through the identity. Differentiation with respect to the parameters x, y, z in the identity provides the generators 



1(7 x 



^M(x,y,z) 



—M(x,y,z) 
dy 

^-M(x,y,z) 
dz 



x—y—z—0 



x=y=z=0 



x—y—z—O 



(10) 

(11) 

(12) 



where M(x, y, z) is defined by Eq. Hence the Lie algebra su(l, 1) is the real vector space spanned by the matrices 
io~ x , i<jy, and a z . By defining k x = i^f-, k y = i-^r, k z = and k± = k x ± ik y , we obtain the standard commutation 
relations 



[k+,k-] = -2k z 
[k z ,k±] = ±k ± 



(13) 



By definition, an operator representation of the algebra su(l, 1) is given by the assignment of three operators K Xl K y 
and K z that satisfy the above commutation relations with K± = K x ± iK y . From such relations, it follows that in 
any representation of stt(l, 1) the Casimir operator 



K-K = Kl-Kl-Kl 



(14) 



commutes with the whole algebra spanned by K x , K y , K z 



B. The exponential map 

A way of writing the group elements in any representation in terms of the Lie algebra generators is through the 
exponential map. The exponential map M = e lm is the map that associates an element m £ su(l,l) of the Lie 
algebra with an element M £ SU(1, 1) of the group. In order to discuss the exponential map, it is suitable to write 
the elements of the algebra as m = \ n • k, where \ £ R, ft ■ k = n z k z — n x k x — n y k y and ft £ R 3 is a normalized 
vector. In this context normalized means that the product ft ■ ft = n 2 z — n x — n y can assume only the values +1, — 1, 
and 0. Then, the exponentiation of the element m £ su(l, 1) is easily performed by using the relation 

(n.fc) 2 = ^l, (15) 

which follows directly from the properties of Pauli matrices. In the following, we analyze the three cases ft ■ ft = ±1, 
separately. 

Case 1: ft ■ ft = +1. The exponentiation gives 

M+ = e lxrri = cos (|) t+i sin (|) 2ft- k . (16) 

Notice that, for any fixed direction n, we have a one-parameter subgroup, which is compact and isomorphic to J7(l). 

The group elements of the form Jl6| ) form a region fi + C SU(1,1), which contains ±1 and all the matrices 
M £ SU(1, 1) such that |Tr[M]| < 2. 
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Case 2: n • n = — 1. Exponentiating the generator ft ■ k we obtain: 

M_ = e lx ™ £ = cosh l+isinh(!) 2n-fc. (17) 

In this case, for a fixed direction n we have a one-parameter subgroup, which is not compact and is isomorphic to R. 
The elements M_ form a region f2_ C SU(1, 1), which contains the identity and all the matrices M £ SU(1, 1) such 
that Tr[M] > 2. 

Case 3: n ■ n = 0. In this case, the exponentiation gives 

M = e 4X "' £ = 1 + i X n ■ k . (18) 

The elements Mq form a region Qq C SU(1, 1), which contains all matrices M G SU(1, 1) such that Tr[M] = 2. The 
region f2o is a two-dimensional surface, and therefore, it has zero volume. 

We want to stress that the exponential map does not cover the whole group SU(1,1). The region of SU(1,1) 
covered by the exponential map is fl = fl + U fi_ U Oqj and contains matrices M £ SU(1, 1) such that Tr[M] > —2. 
However, according to the parametrization (||), the trace of a matrix M £ §U(1, 1) is Tr[M] = 2t, t £ R. Therefore, 
the group §U(1, 1) contains also elements with trace Tr[M] < —2, that cannot be obtained with the exponential 
map. Nevertheless, any matrix M £ SU(1, 1) with Tr[A/] < —2 can be written as M = — M_, for some M_ £ f2_, 
and any matrix M £ SU(1,1) with Tr[M] = —2 can be written as M — —Mo for some Mo £ Qq. Defining 
-Q_ = {-M- | M_ £ and -fl = {-M \ M £ fl } we have 

§U(1,1) =flU-0_ U-O ■ (19) 

Notice that, since fio and — fio have zero measure, any group integral can be written as the sum of only three 
contributions, coming from f2_|_,0_, and — f2_, respectively. 

Even though the exponential map does not cover the whole group SU(1,1), any group element M(6,<p,ip) — 
parametrized as in Eq. (|^) — can be written as a product of exponentials, for example as 

M(6, <j), %j>) = e t k +- ik - e 2i ^ f = -iOe-^-^ . (20) 

The relation ( |20| ) is particularly useful, since it allows to construct from any representation of the Lie algebra 
su(l, 1) a representation of the group SU(1, 1). In particular, for the physical realizations of the group SU(1, 1), where 
the generators k x ,k y , k z are represented by Hermitian operators K x , K y , K z in an infinite dimensional Hilbert space, 



relation (20) provides the unitary representation 

Ubm = e^ K +^ K - e 2 ^ K > £ - -i0e~^~^ . (21) 



C. Baker-Campbell-HausdorfF formula 

The exponential with in Eq. ( ^l|) can be further decomposed according to the Backer-Campbell-Hausdorff 

(BCH) formula. The BCH formula is the fundamental relation, holding for any representation of the algebra su(l, 1), 
given by § 

e iK + -iK- = exp ^r ta " h l«l K + ( -J— ) exp"* tanh l*l K - Ve G C . (22) 

\cosh\Qj 

This formula can be simply proved by verifying it in the case of the two- by- two matrices k+, fc_, k z £ su(l, 1). 
A version of the BCH formula in "antinormal order" is given by the relation 

e iK + -m- = exp-Tf ta»h|e| X_ ( )2 K 2 exp ^ta„h|?| K + \f£ £ <C , (23) 

which follows from (^2|) with the change of representation K' + = —K-, K'_ = —K +1 K' z = —K z . 



III. ELEMENTARY GATES 

The parametrization (|8|) makes evident that any element of SU(1, 1) can be obtained as a product of exponentials 
of the generators k z and k x . In fact, Eq. (^|) is equivalent to the decomposition 

M(0,0,VO = e^-W* e~ 2lk * e^^' (24) 



As a consequence, we have the following approximation theorem: 
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Theorem 1 (Approximation of group elements) Any element of M € SU(1, 1) can be approximated with arbi- 
trary precision with a finite product involving only three elements G\, G 2 , G3 £ SU(1, 1). A possible choice is 

Gi = e 6ia % G 2 = e- 9 ^, G 3 = e i<feff * , (25) 

with 6 1 ,9 2 > 0, 9i/9 2 <£ Q, and 3 /2tt £ Q. 

Proof. Due to decomposition ( p^ ) , it is enough to show that all elements of the form e 1 ^" 7 ' and of the form e 9<Tx can 
be approximated with a product of Gi,G 2 ,G3. First, any point of the circle C = R mod 2tt can be approximated 
by a multiple of an angle ^3, provided that 03 is not rational with 1-k. Approximating <j> as </> « N 3 cj) 3 ,N 3 6 N 
corresponds to approximating the exponential e 1 ^ * as G^ 3 . In the same way, any point of the circle C = K mod #1 
can be approximated by a multiple of — 9 2 , provided that 9 2 is not rational with 9%. Since any real number S£l can 
be written as 9 — M9\ + 9 mod 81, by approximating 9 mod 9\ ~ mod 8\, we obtain w N\9x — N 2 6 2 , for 

some iVi 6 N. This corresponds to approximating the exponential e 6 * * as G^G^ 2 . 

The previous theorem is particularly important in consideration of physical realizations, where the group SU(1, 1) 
acts unitarily on an infinite dimensional Hilbert space. In this case, the previous result shows that any unitary 
transformation representing an element of SU(1,1) can be arbitrarily approximated by a finite circuit made only 
of three elementary gates. However, if we thoroughly define a parameter for the rating of the approximation we 
find that the accuracy is arbitrarily small, and this fact is due to unboundedness of the generators for the unitary 
representations of SU(1, 1) of physical interest. In particular, we are interested in the two representations in which 

K z = -(ata + tfb + 1), K+ = a i b i , K- = ab, 

\ ! ! (26) 

K z = -(ota + 1/2), K + = -{a))\ K_ = -a 2 . 

The parameter for the approximation rating is the accuracy e _1 , with 

e= sup KUi-UtMl (27) 
11-011=1 

where U± is the target element, U 2 is the product of elementary gates that approximates U\. However, since we 
are considering infinite dimensional representations, the difference U\ — U 2 has eigenvalues arbitrarily near 2. The 
supremum is then always 2, and in order to find some approximation criterion one has to impose some constraint on 
the states that we are considering. For example, we will impose that the average and second moment of the photon 
number distribution are finite, which are reasonable physical assumptions. Suppose now that we have a sufficiently 
long sequence of elementary gates, in such a way that, using the decomposition of Eq. (0), U\ — e~ %aK " e~^ Kx e~ % ' iK ' s 
and U 2 — e~ l ( a + s °-) K z e -<-{P+ s i3) K x e - l {'i+ s ~t) K z j anc i orU y nrs t order terms in S x are relevant, thanks to the constraint 

on states. After some algebra and exploiting Eq. ( |63| ) one can verify that the supremum of (if>\21 — U\lJ 2 — U\U\\^>) 
is almost equal to the supremum of (^|A|^>), where 

A = f (S 2 a + 5* + 2 co S h(3S Q 5, - 5l)K 2 z 

\(^-^-^-2 cosh/3<U 7 )iC ' " S ' 

depending on the sign of b\ + S 2 + 2 cosh/3<5 Q (5 7 — 6 p. This equation implies two facts. First of all, we can easily 
verify that the physical constraint on states is necessary in order to guarantee boundedness of e. However, it is not 
sufficient because of the presence of cosh/3 in the expression. This is due to non compactness of the group, which 
implies that even in the defining representation the approximation is worse as one goes further along the direction of 
a non compact parameter. This fact fatally flaws any analogy to the Solovay-Kitaev theorem for the qubit case, and 
in order to have a similar result one must also restrict the set of unitaries that he wants to approximate. Otherwise, a 
power law for the number of gates as a function of e _1 can be searched which contains an explicit dependence also on 
the parameter (3. Suppose that we have |<5„ + S? ( + 2cosh/3<5 a <5 7 — <5l| = f(N), where N is the number of elementary 
gates needed to approximate the target group element with in the defining representation. Then, for A oc K 2 in 
Eq. d28|), we have e = f(N)((E 2 ) + 2X(E) + A 2 ), where E is the total number of photons and A = \ for single mode 
representation and A = \ for two modes. The function / is clearly non increasing, and supposing that it is strictly 
monotonic, it can be inverted, obtaining 



N=r i m+tm+v ■ (29) 
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IV. UNITARY REPRESENTATIONS OF SU(1,1) 

Given a representation of the su(l, 1) algebra where the generators K x , K y , K z are Hermitian operators acting in an 
infinite dimensional Hilbert space TL, we consider the unitary representation U$ ^ of the group SU(1, 1) defined by Eq. 
(pl|). In general, such a representation is reducible, and it can be decomposed into unitary irreducible representations 
(UIRs). 

A UIR Ug^^ is called square- summable if there is a nonzero vector |t>) G TL such that 

6v(6,<f>,^) U 0>M \v)\ 2 < oo , (30) 

/SU(1,1) 

where dv is the invariant measure defined in Eq. @. Moreover, since the group SU(1, 1) is unimodular, if the above 
integral converges for one vector \v) ^ 0, then it converges for any vector in TL JTof . 

Square-summable representations enjoy the important property expressed by the following: 

Theorem 2 (Formula for the group average) If the irreducible representation Ug^^ is square-summable, then 
for any operator A £ B(TL) the following relation holds 

dv{9, 4>, i>) Ug, M A = Tr[A] ~ . (31) 

SU(1,1) a 

Here 1 is the identity in TL, and d is the formal dimension, defined by 

d=([ dv(9,4>,i>)\(v\Uo,M\v)A , (32) 
\ysu(i,i) J 

where \v) is any normalized vector \v) € TL, (v\v) = 1. 

The formula for the group average is fundamental in the contexts of quantum estimation and tomography, since it 
allows to construct resolutions of the identity via a group integral. In the context of quantum estimation, Eq. (|3l| ) 
ensures that the operators 

P(6,4>,il>) = U e ,Mtut >Mt (33) 

where £ is any operator satisfying £ > 0, Tr[£] = d, provide a positive operator valued measure (POVM) for the joint 
estimation of the three parameters 9,(f>,ip. In fact, such operators satisfy the normalization condition 

dv(6, 0, VO P{9, 0, = 1 , (34) 

SU(1,1) 

which guarantees that the total probability of all possible outcomes is one. In particular, if £ = d \v){v\ for some state 
\v), the above formula gives the completeness of the set of SU(1, 1) coherent states O] 

\ve,<t>,i>) = Ue,,},,-,!, \v) . (35) 



A. Examples 

1. Single mode squeezing 

The representation of the su(l, 1) algebra, given by 

a t2 a 2 1 / t 1\ 

K + = - K. = - K n = -(cta+-) (36) 

is reducible in the Hilbert space H of a single harmonic oscillator. In fact, the subspaces Tt cvcn = Span{|2n) | n € N} 
and 7i dd = Span{|2n +1) | n € N}, defined in terms of the Fock basis |n) = -^aT n |0), are invariant under the 
application of K x , K y , K z . The unitary representation of SU(1, 1) defined by Eq. acts irreducibly in the subspaces 
7i e ven and TLodd- There is a substantial difference between the two UIRs acting in 7i dd and Ti cvcn , in fact, the first is 
square-summable, with formal dimension <i dd = 1/(4tt 2 ), while the latter is not. A surprising consequence of the non 
square-summability in H. cvcn is that the squeezed states |£) = e^ K+ ~^ K ~ |0) — which are the coherent states of SU(1, 1) 
commonly considered in quantum optics — do not provide a resolution of the identity. 
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2. Two modes squeezing 



The representation of the Lie algebra su(l, 1) given by the operators 

K+=o)b\ K_=ab K z = i (0*0 + 6*6 + 1) (37) 

is reducible in the Hilbert space TL a ® TLb of two harmonic oscillators. It is indeed immediate to see that, for any 
5 G Z, the subspaces TLs = Span{|m)|n) | m, n S N, m — n = <5 } are invariant under application of the operators 
K + ,K^ 1 K Z . The unitary representation of SU(1, 1) given by Eq. ( fjl| ) is irreducible in each subspace TLs- The two 
UIRs acting in TLs and TL^s are unitarily equivalent, while for different values of \5\, one has inequivalent UIRs. All 
UIRs in the two modes realization are square-summable, with the only exception of the case 6 = 0. 



V. SU(1, 1) TOMOGRAPHY 



A. Reconstruction formula for square-summable representations 



Let us now consider the group as a tool for quantum tomography. In order to do that, it is useful to consider the 
set of operators on a Hilbert space Ti as a Hilbert space itself, isomorphic to Ti ® Tt, and to look for spanning sets 
in this Hilbert space. An immediate and handy way of defining the isomorphism between operators and bipartite 
vectors is through the definition 

\A)) = ^2(m\A\n) \m) ® \n) (38) 

where A is an operator on TI, and \n) are elements of a fixed basis for Ti. This definition implies the following useful 
identities 



A®B\C)) = \ACB T )) (39) 
((A\B)) = Tv[A^B] , (40) 

where X T denotes the transpose of X in the basis |n). 

For a square-summable UIR Ug^^, we can obtain a resolution of the identity by simply exploiting Eqs. ([n]) and 
(|it), namely 



l®l = d/ dv(6,<i>,ip) \U e , M ))((U e ,^\ , (41) 

which shows that the unitaries Ug^^ form a spanning set for the space of operators. 

Tomographing the state p is equivalent to reconstructing the ensemble average Tr[pyl] = ((p\A)) of any operator A 
on the state p. This can be done using the reconstruction formula 

Tr[pA] = d [ M9, 0, V) Tr[pUg^} Tr[^ , , (42) 

JSU(l.l) ' ' 

which directly follows by inserting the resolution of the identity ( [4l| ) into the product Tr[pA] = ((p\A)). In a real 
tomographic scheme, the traces Tr[pUg t< f,^] have to be evaluated by experimental data, and subsequently averaged 
with the processing function /^(0,<^, tp) — Tr[?7g ^ ^A] in order to obtain the expectation value Tr[pyl]. A feasible 
scheme for evaluating the traces Tt[pUg^^\ by experimental data is discussed in Subsection [~~ 



B. Non square-summable representations: regularization 



The representations of SU(1, 1) that ar e comm on i n quan tum optics are single mode and two-modes Schwinger 
representations, analyzed in Paragraphs. [V A 1 and IV A 2 , respectively. In particular, the irreducible subspaces 
Heven in the single-mode case, and TLq in the two-modes case are particularly interesting, since coherent states with 
even photon number in the single-mode case (or, alternatively, zero difference of photon numbers in the two-modes 
case) are experimentally achievable by simple vacuum squeezing. 
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The problem now is that, the single-mode representation in Tl cvcni and the two- modes representation in TIq are not 
square-summable, therefore the group integral in Eq. ( |30|) diverges. In order to circumvent this problem, we address 
here the a technique that consists in modifying the invariant measure d^(0, </>, iji) by a regularization factor g(9, 4>, if;), 
which is positive almost everywhere. The modification of the measure makes it non invariant, and consequently the 
group average identities become similar to those of non unimodular groups, where there is no invariant measure. 

Using the regularization factor, instead of the resolution of the identity ([4l|), we have a positive invertible operator 

F= I di/(M,^) 0(M,1/O \Ue, M ))((U e>M \ , (43) 

«/S0(l,l) 

and the ensemble average of any operator A can be obtained by writing Tr[jaA] = ((p\FF~ 1 \A)) . In this way, we can 
provide a regularized reconstruction formula 

Tr[pA}= [ di/(M > V0ff(M > V0/A(M,V0^[tfM,M ( 44 ) 
involving the processing function 

fA(8,M) = ({Ue, M \ F ~ 1 \ A h (45) 

instead of ((Ug^^\A)}. Notice that the identity Eq. ( f44| ) can be used also in the square-summable case with g(6, <f>, ip) = 
1. 



C. A relevant example 

Here we consider in detail the case of the UIRs with even photon number (single- mode case), and with 5 = 
(two-modes case) representations, providing an example of the general method discussed above. We will start from 
the following integral 

S(m,n;m',n')= [ di/(0, 0, iP)(2m\U e .^\2n)(2n'\ul . J2m'), (46) 
Jsv(i.i) ' ' 

where |2m) denotes both the even eigenstate of a^a or the zero-difference eigenstate \m)\m) of a) a + b^b. This 
can be evaluated by exploiting Eqs. (^l|) and thus obtaining the following expressions for the matrix element 

(2m|£/0»|2n) 

(2m\Ug^\2n) = e <*(2r.+K) e <W-*X»-f») £ e K (p)(-*tanh0) 2p+m -" (— —J , (47) 

p=0 \COS J 

where k = 1/2 for even single mode states and k = 1 for d = two modes states, and 



V2n!2m! 1 

K = — 



c K (p) = { pl(j) + m ~ " )!(2n ~ 2p)!2 2 f+ m ~" 2 (4g) 

n\m\ 

K = 1. 



p!(p + m — n)!(n — p)! 2 
By exploiting the integral in ip and then in </> we obtain 



where 



S(m,n- : m ,n) = Air 2 5 m ,m'5n,n' ^ c K (p)c K (p')(-l) p+p I K (m, n,p,p'), (49) 



I K (m,n,p,p>) =/ d0sinhflcoshfl \ - , (50) 



(cosh6») 4n - 2 P- 2 P'+ 2K ' 
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and for 2n = p + p' this is clearly divergent. Moreover, for m < n and p = p' = the integral diverges because of the 
singularity in 9 = 0. If we introduce the regularization factor g(0,<p,tp) = g(9) = e ( coshS yi — , by exploiting the same 
calculations we get the same result with I K (m,n,p,p') substituted by 

I K . g (m,n,p,p') = f dxx 2n - p - p ' +K (l~x) p+p ' +m - n e'^, (51) 
Jo 

which is derived from Eq. ( f30| ) by the change of variable (l/cosh#) 2 — > x, and which is finite. As a consequence 

oo 

dv{d,ct>,^)g{d)\U e ^)){{U B>M \= F^ n \2m)(2m\®\2n){2n\ 1 

m .n=0 

n 

^1=4tt 2 j2 {-iy+P'c K (p)c K (p')I K , g {m,n,p,p') (52) 

p,p'— 

1 - r\ p 2 

' s x 2n (l-x) m - ,l e- — , 



l 

d.c 



(i 



^(-l) p c K (p) 



and since the coefficients c K (p) are non null, clearly < Fmli < oo. This implies that the operator F in Eq. (|42 
is actually invertible, and we can safely use the processing function in Eq. fliq ) for tomographic reconstruction of 
the operator A. Notice that this formula is very close to the group-average formula for non unimodular groups, 
where the Duflo-Moore operator C is involved. The group average identity in that case is similar to Eq. fljffi ) with 

Fm,n — (C^C)m,m- 

D. How to make tomography experimentally 

In quantum tomography, the expectation value (A) = Tr[pA] of any observable is reconstructed by exploiting 
integral (pi). Moreover, in order to have a feasible tomography, it is essential to devise a method to evaluate the 
traces Tr '[pUe^^] from experimental data. To do this, it is useful to break the integral over SU(1, 1) into the sum of 



the contributions coming from the regions f2 + ,5!_ and — fi_, introduced in Par. (IIB). It is not difficult to see that 



the regions Q- and — f2_ give the same contribution to the tomographic integrals, whence we have 

Tr[pA] = f di/(M,^)$(M,V0/A(M,V0'&[tf9,*,*P] ( 53 ) 
Jn+ 

+2 f dv(6, 0, V) 9(0, <P, WaV, <P, i>)Tt[Uo, MP ] . 



By definition, any element in J7 + (^-) can be obtained by the exponential map as e txn ' K for some x an( i ^ with 
ft ■ n — +1 (—1). In addiction, it is possible to show that any exponential e %xn ' K can be written as 

j X n-R = / V(n)l e*** V(n) , ft ■ ft = +1 , , 

\ W(ny e ixK * W{n) , n-n = -l { ' 

where V(n) and W(n) are suitable unitaries in the group representation. A detailed proof of this result is given in the 
Appendix. Thanks to this observation, the trace Tr[pe txn K } can be evaluated by performing a unitary transformation 
on the state p (either V(n) or W(n)), and subsequently by measuring one of the observables K z and K x . 

Finally we observe that, since a real experiment produces only a finite array of data, the integral ( |5^ ) has to be 
be approximated by a statistical average over the experimental results obtained by measuring a large number N of 
identically prepared systems. This introduces the need of a randomization in the experimental setup, that produces the 
unitaries V(n),W(n) according to some probability distribution. Notice that the most natural choice, that would be 
to take dp(n) as the measure over the space of directions ft induced by the invariant measure du(9, tfi, ip) is not possible, 
since such a measure cannot be normalized (the space of directions is noncompact). The form of Eq. ( |53|) suggests 
then to take as a measure dis(9, </>, ip)g(9, </>, ip), and in the example we considered this actually works. However, it may 
happen that regularizing the integral in Eq. (^) is not sufficient for regularizing also the group measure. In this case 
it is convenient to modify g(9, 4>, tp) in such a way that both the measure itself and the group integrals converge. This 
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implies in particular that the choice g(0, cf>, iji) = 1 for square-summable representations has to be changed. Finally, 
the ensemble average (A) can be then be approximated by the expression 

1 N 

< A > ^ E /^i.^.^W^.^,^)]. (55) 
i=i 

where are the randomly extracted parameters. Notice that the expression on r.h.s. in Eq. (|55j) reasonably 

converges to l.h.s. if the variance of the processing function is finite, namely if fA(0,(j),ip) is square summable. By 
Eqs. ( |43| ) and (jij]) this condition is equivalent to 

/ dKM,#?(M,^)IMM,VO| 2 = ((AlF-^A)) < co. (56) 

•/SU(1,1) 



VI. CONCLUSIONS 



This paper collects a large number of useful results about the group SU(1, 1) that are dispersed in the literature, and 
also contains some novel applications regarding the use of SU(1, 1) for quantum computation and tomography with 
nonlinear optics. The main issues we addressed here are i) the approximation of SU(1, 1) gates in the quantum optical 
representations and ii) the tomographic state reconstruction exploiting group theoretical methods. As regards the 
first topic, we gave an approximability theorem and discussed the limits under which it holds. The theorem provides 
a useful result in the search for an elementary set of gates that can be used to universally approximate any §U(1, 1) 
gate with arbitrary accuracy. To complete the analogy with the Solovay-Kitaev theorem for qubit gates, the power 
law of the number of elementary gates as a function of the accuracy should be evaluated, and due to non compactness 
we expect that the law would depend on the parameters of the target group element. 

In the context of quantum estimation and tomography, we showed a technique for regularization of the group integral 
for the physically relevant representations, that are not square-summable. The core of the regularization technique 
is a modification of the Haar measure over the group, such that the regularized measure is no longer invariant. This 
makes the integrals for tomographic reconstruction convergent, but radically modifies the processing functions. Such 
a regularization technique is very powerful, since it contains a freedom in the choice of the regularization factor, that 
allows for a further optimization of the processing. Moreover, the mentioned scheme can be applied not only to the 
case of SU(1, 1), bus also to any other tomographic setup. 
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VII. APPENDIX 




For a given representation of the su(l,l) algebra, consider the real vector space V spanned by the generators 
K x , K y ,K z . Of course, V is isomorphic to K 3 via the correspondence 

K y «-► (l\ K z «-> . (57) 

The action of the group SU(1, 1) on the space V, given by V 3 m 1— > e %xn ' K m e~ %xn , can be obtained by exponen- 
tiating the adjoint action on the algebra, namely 

e ixH - 6 m e - ixH g = e^'AdW m , (58) 

where Ad(ifj) is defined by Ad(Ki)Kj = [Ki,Kj]. Moreover, using the commutation relations of su(l, 1) it is 
immediate to find that 

/o o\ /00A /0 -i o\ 

Ad(K x ) = -i\ Ad{K y ) = Ad(iQ = \i . (59) 
\0 -i / i 00 \0 0/ 
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Therefore, we obtain that a generic element of SU(1,1) — parametrized as M(9,<p,i/j) — e 1 ^ ^ kz e 2lk * e lly ' t ' + ^ kz as 
in Eq. (g4|) — is represented in the space V by the matrix 

R(9, cb, j>) = e <(*-*)Ad(Jf.) e -WAd(Jf.) e i(<f>+4)Ad(K z ) (60) 

whose explicit expression is rather lengthy, but easily computable by exponentiating the matrices in Eq. (p9[). 

It is not difficult to see that the matrices R(9, cj>, -0) given by Eq. (6(J) form a subgroup of the group §©(2,1), 
namely they all have unit determinant and preserve the form v ■ w = v z w z — v x w x — v y w y . More precisely, the matrices 
R(9, <f>, ip) coincide with the group SO + (2, 1), which contains all matrices R € SO(2, 1) such that R33 > 1. Incidentally, 
we notice that the correspondence SU(1, 1) — * SO + (2, 1) is not one-to-one, since both ±1 £ SU(L 1) are mapped into 
the identity in SO + (2, 1). One has indeed the group homeomorphism SO + (2, 1) ~ SU(1, 1)/Z2[jl4|, which is exactly 
the same relation occurring between the groups §U(2) and SO(3), namely SO(3) ~ SU(2)/Za. 

Similarly to the case of SO(3), where any spatial direction ft can be conjugated with the direction of the z— axis 
by a suitable rotation, in the case of SO + (2, 1) any direction ft with ft ■ ft = +1 can be conjugated with the z— axis, 
and any direction with ft • ft = —1 can be conjugated with the a;— axis. For example, the matrix R(9, <f>, m Eq. ( |60| ) 
transforms the direction of the z— axis as 

'-sinh(20) sinO-i/^ 

- sinh(26>i !cos(0 - tp) | , (61) 
cosh(26») 

and it is clear that here ft can be any direction with ft ■ ft = 1 (modulo an overall phase factor). Therefore we have, 
for any ft with ft ■ ft = +1 



n-K = K z Ue i<Pii , , (62) 



for suitable 9, 0, ip. In conclusion, 



'6,<t>,-<t> 

The same reasoning holds for any direction with ft ■ ft, — — 1 



xn ' K = Ul M e %xK * Uo,w ■ (63) 
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